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Abstract 



We study the effect of external stochastic modulation on a system with 
0(2) symmetry that exhibits a Hopf or oscillatory instability in the absence 
of modulation. The study includes a random component in both the con- 
trol parameter of the bifurcation and in the modulation amplitude. Stabil- 
ity boundaries are computed by either solving the stationary Fokker-Planck 
equation in the vicinity of the center manifold of the underlying deterministic 
system whenever possible, or by direct numerical solution otherwise. If the 
modulation amplitude has a stochastic component, the primary bifurcation is 
always to standing waves at a value of the control parameter that depends on 
the intensity of the fluctuations. More precisely, and to contrast our results 
with the case of a deterministic periodic forcing, the onset of instability in 
the standing wave regime is shifted from its deterministic location, and the 
region of primary bifurcation to traveling waves disappears yielding instead 
standing waves at negative values of the control parameter. 
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I. INTRODUCTION 



At a Hopf bifurcation in a periodically modulated system, the trivial state loses stability 
to either traveling or standing waves above onset depending on the amplitude of the mod- 
ulation b. For sufficiently small modulation amplitudes, traveling waves appear at a fixed 
value of the control parameter, or, independent of b. The threshold for standing waves, 
however, is a decreasing function of b. We discuss in this paper how the existence of a 
stochastic component in both cir and b affects the nature of the bifurcation, as well as the 
stability boundaries of the trivial state. The calculations presented here are not specific to a 
particular system, but rather are based on the normal form equations appropriate for a Hopf 
bifurcation in a system with 0(2) symmetry when driven by a periodic force of frequency 
about twice the Hopf frequency of the unperturbed system. 

Detailed studies of Hopf bifurcations have been given for a large number of systems [|TJ. 
We mention, for example, the transition from straight rolls to Busse oscillations observed 
in Rayleigh-Benard convection. This instability occurs in fluids of low Prandtl number (P r ) 
and at sufficiently large values of the Rayleigh number (R). For instance, it is observed in 
air (P r = 0.71) when R reaches a value close to 6000. The instability manifests itself as a 
periodic transverse distortion of the rolls that propagates along their axes. Recently, Clever 
and co-workers |2j have studied the influence of a periodic modulation of the gravitational 
field on the instability. To that end, they solved numerically the time-dependent nonlinear 
equations for three-dimensional convection using values of R above the onset of oscillatory 
convection. They varied the amplitude and frequency of the modulation, with the latter 
always set to a multiple of the fundamental (unmodulated) frequency. The main result 
of their study is that for moderate values of the modulation amplitude, a transition from 
traveling to standing waves takes place, with the system's frequency response being either 
synchronous or sub-harmonic. The authors also found that, as the amplitude is further 
increased, this frequency locking behavior disappears and a time-dependent aperiodic regime 
sets in. 
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The onset of convection in binary fluids also occurs through an oscillatory instability 
when the separation ratio is negative, i.e., when the temperature field is destabilizing whereas 
the composition gradient is stabilizing Given the large difference in time scales between 
energy and mass diffusion, the process is also known as double-diffusive convection. This type 
of instability is commonly observed in directional solidification experiments when a crystal 
which is being grown upwards rejects a heavier solute. The effect of a periodic modulation of 
the gravitational field has been addressed theoretically by Saunders et al. [[| for a laterally 
unbounded fluid layer and stress free boundary conditions at the top and bottom of the 
layer. In the region of parameters in which the bifurcation of the unmodulated system 
is oscillatory, they find that below onset of the unmodulated system there exist regions of 
instability to stationary convection (either subharmonic or synchronous with the modulation 
frequency) for sufficiently large values of the modulation amplitude. For conditions above 
onset of the unmodulated system, the bifurcation is to traveling waves for arbitrarily small 
amplitudes of the modulation. These findings are in agreement with the general bifurcation 
diagram for a system of 0(2) symmetry that will be discussed below. 

Experiments on double-diffusive convection are often performed between two conducting 
plates which are also made permeable so as to control the solute concentration at the top 
and bottom of the system. As in the Rayleigh-Benard case, the primary instability leads to a 
pattern of rolls, but that now propagate through the system. Experiments performed in long 
narrow cells of annular or rectangular geometry have shown a traveling wave pattern that 
is either uniform in space or confined to a small region of the cell [|J ||. The influence of a 
periodic modulation of the temperature gradient has been investigated by Rehberg et al. 
who studied convection in a water-ethanol mixture in a small rectangular cell, starting from 
a uniform traveling wave pattern. As was the case in the gravitationally modulated fluid 
layer, they observed the emergence of a standing wave structure as a periodic modulation 
of sufficient amplitude was added. 

In the same paper, these authors report results from a much more elaborate study of the 
onset of electro-hydrodynamic convection in the nematic liquid crystal Merck Phase V. In 
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this system, an electrostatic potential difference applied across the experimental cell plays 
a role similar to that of the temperature gradient in thermally induced convection. In order 
to suppress charge injection processes at the electrodes, the applied voltage is alternating at 
a frequency to. As its r.m.s. value is increased, the motionless state loses stability to a roll 
pattern, the properties of which depend on the driving frequency u. For the Merck phase V 
system, steady Williams rolls emerging at low frequencies give way to spatially homogeneous 
traveling waves as the driving frequency is increased. The authors studied the stability of 
these traveling waves against a small periodic modulation of the voltage. This perturbation, 
superimposed on top of the basic ac driving, had a small frequency uj m in resonance with 
that of the traveling waves (i.e., uj m ~ 2u TW « w, with uitw the frequency of the traveling 
waves in the unmodulated state). As in the two cases described above, traveling waves were 
found to lose stability with respect to standing waves as the amplitude of the modulation 
was gradually increased. A corresponding shift in the threshold was observed, with the 
convecting state appearing at smaller values of the control parameter (i.e., the r.m.s voltage 
of the ac source). 

A general description of a Hopf bifurcation in a periodically modulated system has been 
given by Riecke, Crawford and Knobloch J7J. Their analysis, which is briefly reviewed in 
Section |T[ involves two complex amplitude equations governing left and right traveling 
waves emerging at a Hopf bifurcation. The periodic modulation, which is assumed small, 
provides a linear coupling between the two, and leads to the excitation of standing waves 
under certain conditions. Different branches of the bifurcation diagram mark the onset 
of standing or traveling waves, and they join at a codimension-2 bifurcation point which 
has been observed in the electro-hydrodynamic convection experiments of Rehberg et al. 
The model also predicts a number of secondary instabilities which have yet to be observed 
experimentally. 

The purpose of this paper is to extend the results summarized above to cases in which 
either the control parameter or the amplitude of the modulation fluctuate randomly. Al- 
though this is a quite general question, we are specially motivated by experiments conducted 
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in a microgravity environment 0|§. There, the effective gravitational field is known to fluc- 
tuate in time with the amplitude of the fluctuation being two or three orders of magnitude 
larger than the residual steady gravitational field [H| . The frequency spectrum of the resid- 
ual acceleration field, or g-jitter, typically comprises periodic components and a white noise 
background |TT]]. The physical origin of these disturbances lies in the many mechanical 



processes that take place onboard spacecraft, and their coupling to mechanical modes of 
the structure. A recent analysis of actual acceleration data taken during a Space Shuttle 
flight has shown the existence of several periodic components of frequencies in the range of 
a few Hz, and amplitudes of the order of lO -3 ^, where $e is the intensity of the Earth's 
gravitational field. There appears to be also a white noise background with approximatively 
gaussian statistics. We attempt to present here a general framework within which to ana- 
lyze the effects of such a residual field on an oscillatory instability, and therefore to provide 
the basis for future studies of specific systems. Two areas of concern include the appear- 
ance of undesired instabilities of some base state caused by g-jitter, and the modification 
in character and location of onset of a given instability because of the random component 
of the effective gravitational field. The specific cases of directional solidification and double 
diffusive convection under reduced gravity conditions have been reviewed in W2\. 



In the classical deterministic system is said to undergo a bifurcation when its 

long time behavior changes qualitatively as some control parameter is continuously varied. 
Mathematically, this change corresponds to an exchange of stability between different so- 
lutions to the system's governing equation(s). The nature of the bifurcation depends on 
that of the solutions it involves: the saddle-node, transcritical and pitchfork bifurcations, 
for instance, all involve two fixed-point solutions, while the Hopf bifurcation has both a 
fixed-point and a limit cycle. Each one of them has an associated set of equations, known 
as its normal form, to which any specific example transforms in a small region around its 
bifurcation point. The dimension of this set is equal to the smallest number of equations 
that can still give rise to the bifurcation (one equation in the first three examples given 
above, two in the Hopf case). In systems described by a larger number of equations than 



those involved in the normal form, the governing set of equations can be reduced close to 
the bifurcation point. The reduced set defines a surface in the phase space of the original 
equations known as the center manifold. The existence of this surface, which has the same 
dimension as the normal form, therefore leads to a simplified formulation of the problem, 
with an underlying separation of time scales in the evolution of variables on, and orthogonal 
to, the center manifold. 

The effect of random fluctuations, both of internal and external origin, on bifurcations 
has been studied in considerable detail []13|-[16[|. Internal fluctuations, typically of thermal 
origin, enter the governing equations linearly or "additively" , scale with the inverse of the 
system's size, and lead in general to so-called imperfect bifurcations: the bifurcation point 
is smeared into a small region of size proportional to the intensity of the fluctuations. On 
the contrary, externally induced fluctuations (e.g., random changes in the externally set 
control parameter for the bifurcation) typically enter the governing equations nonlinearly or 
"multiplicatively" , do not satisfy any a priori scaling with the size of the system, and the 
bifurcation point may remain sharp, although its position can depend on the intensity of 
the noise. For instance, Graham |TB| has shown that a Hopf bifurcation with a fluctuating 
control parameter exhibits a sharp rise from the trivial state, without any shift in the location 
of threshold, and a small decrease in the average value of the amplitude of the new stable 
state compared with the deterministic case. As will be shown below for white gaussian 
noise, these results also hold when a periodic modulation of the control parameter is added 
to the system. If the intensity of this modulation is also allowed to fluctuate however, the 
system's response is much more complex. It leads, among others, to shifts in threshold and 
to excitation of standing waves in a region of parameters in which they were previously 
absent. 

The remainder of this paper is structured as follows: Section [TI] briefly reviews known 
results on the effect of a resonant modulation on a Hopf bifurcation in a system with 0(2) 
symmetry. The stochastic extension of the analysis is given in Sections |TJ and |TV[ The 
former discusses the case of a stochastic component in the control parameter an while the 
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latter addresses the case of a random forcing b. 



II. HOPF BIFURCATION UNDER PERIODIC MODULATION 



The results presented in sections |TJ and [TV] extend the work of Riecke and co-workers 



on Hopf bifurcations in periodically driven systems [^]. For completeness, we present a brief 
overview of their work here. Close to a Hopf bifurcation, two complex amplitude equations 
are needed to describe the slow evolution of the unstable modes. Let \1/ describe the state 
of the system. Then, 

^ = Ul {t)e iqz + u 2 (t)e iqz + c.c. (1) 

The complex amplitudes u\ and u 2 correspond to the two eigenvalues ±iwn associated with 
the bifurcation (with wh the Hopf frequency of the limit cycle), and q is the characteristic 
wavenumber of the emerging structure (inversely proportional to the roll width in convection 
experiments). General equations governing the evolution of u\ and u 2 are obtained by 
imposing their invariance under both spatial translations (T : z — > z + d) and spatial 
reflections (K : z — > —z) (0(2) symmetry). From Eq.(|T]), we have T(u±, u 2 ) = (e iqd ui, e tqd u 2 ) 
and K(ui,u 2 ) = (u^ul). Since the equations for ui and u 2 must remain invariant under 
these transformations, they have the form 

d t Ui = g x u\ + g 2 u 2 

d t u 2 = g* 2 ui + g{u 2 , (2) 

where g\ and g 2 are nonlinear functions of the invariants |iii| 2 + |w2| 2 , u\u 2 and u{u 2 , and of 
the external modulation a. The analysis is further restricted to the strong resonance case 
in which the frequency of the modulation is almost twice the natural frequency wh of the 
system, or a ~ 5 e 2 ™H* ; with b real. Letting u± = t]e tWHt and u 2 = £e~ lWHt in Eq. (g) and 
dropping oscillatory terms, one obtains to cubic order the equations 

d tV = ari + b( + cr/(|r/| 2 + |C| 2 ) + gv\C\ 2 (3) 



and 

d t c = a*c + b v + c*((\v\ 2 + ICI 2 ) + gX\v\ 2 , (4) 

governing respectively the evolution of left (rj) and right (() traveling waves. The real part 
of a (a R ) is the control parameter while its imaginary component (a«) is the detuning of the 
wave from subharmonic resonance. The stability of the trivial state i] = ( = is determined 
by linearizing Eqs. (|3|) and (^) and letting rj = fje and ( = Qe xt . Solutions of that form 
exist provided 



- brj I 

C = - and A = a R ± Jb 2 - of. (5) 

A — Q 



Thus, the system undergoes either a steady bifurcation at clr = ±y& 2 — of (b > aj), or 
a Hopf bifurcation at a# = (b < aj). The point (apt,b) = (0,<2j) delimiting the two 
corresponds to a codimension-two Takens-Bogdanov (TB) bifurcation point. 

In order to study the full non-linear behavior of the system, it is useful to introduce the 
notation r] = xe lipi , ( = ye lip2 , x — <Pi ~ ¥2 and <fi = <pi + <f2- Then, 

d t x = a R x + by cos \ + c R x{x 2 + y 2 ) + g R xy 2 (6) 

d t y = a R y + bx cos x + c R y(x 2 + y 2 ) + g R yx 2 (7) 

d t x = 2di + ni(x 2 + y 2 ) - b sin x(x 2 + y 2 ) / xy , (8) 
with rii = 2ci + gi. The phase angle 4> obeys the decoupled equation 

d t (f> = bsmx(x 2 - y 2 )/xy - gi(x 2 - y 2 ). (9) 
In terms of these new variables, the state of system \l/ now reads 

Equations (|§),(0) and (||) admit two types of stationary solutions: standing and traveling 
waves. For standing waves (SW) x = y and d t <p = 0. In that case, 
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J , =v , = _ Jtf i±[i-^(4+ a ?-y)/^ i 

with M = aitii + aRiiR, N 2 = n R + n 2 and = 2cr + g^. For traveling waves (TW), x ^ y 
and 7^ 0. They correspond to solutions 



•i-- ... 



-a fl [l±(l-4A 2 ) 1/2 ]/2c fl , y r ,, = ar /jr> (12) 



with A 2 = b 2 c 2 R /(a 2 R g R + 4f2 2 ) and = ajCR — a r rii/2. These solutions exist as long as 
A 2 < 1/4 at which point the left and right traveling waves merge to form a standing 
wave. The solid lines in figure [3] delimit the various regions of the stability diagram for the 
parameter set (oj = 2, cr = —1, q = 2, g R = —1, g± = 1). 

In summary, for small modulation amplitudes the system behaves exactly as in the 
unmodulated-modulated case: traveling waves appear at onset, which is located at clr = 0. 
For modulation amplitudes larger than the detuning, standing waves are excited instead, 
the threshold is at or < 0, and is a decreasing function of the modulation amplitude b. 

III. STOCHASTIC MODULATION OF THE CONTROL PARAMETER 

We begin this section with some brief considerations about the study of bifurcations 
in a stochastic system. As already mentioned in the introduction, we do not consider 
fluctuations of internal origin (thermal fluctuations, for example), but rather fluctuations in 
the externally set control parameters. The latter are not necessarily small, typically enter the 
equations nonlinearly or "multiplicatively" , and their effect is not generally a simple smearing 
of the deterministic threshold (the so-called imperfect bifurcation in the case of fluctuations 
of internal origin). Leaving aside the mathematical complexity involved in treating but the 
simplest cases, there remains some discussion in the multiplicative case about such basic 
questions as the proper definition of the threshold, or the degree of generality of the results 
obtained vis a vis the particular details of the model equations or the statistical properties 
of the fluctuating components. As there is no general agreement on these issues, we first 
outline our underlying assumptions here. 
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Our definition of instability or bifurcation point follows the work of Graham |13| , and is 
based on the stationary solution of the Fokker-Planck equation for the system of interest. If 
the system bifurcates from the trivial state, the solution below threshold is a delta function 
centered at zero. The onset point corresponds to the value of the control parameter at which 
additional stationary solutions of the Fokker-Planck equation appear with some nonzero 
moments. Other non-normalizable stationary solutions that may appear below this onset 
are not considered. Second, the dimension of our starting set of equations (||) and (|J) 
is larger than that of the unstable manifold of the deterministic case. We have adopted 
a center manifold reduction procedure in the stochastic case which is analogous to the 



one proposed by Knobloch and Wiesenfeld [17]. The stationary probability distribution 



function is assumed to factor into a contribution that depends only on the slow variables, 
and another that confines the evolution of the system to a small region around the center 
manifold of the underlying deterministic system. Numerical evidence is presented supporting 
such a factorization. We note that, unless fluctuations in the direction normal to the center 
manifold can be completely neglected, such a procedure is not equivalent to adiabatically 
eliminating the fast variables directly from the original model equations, a procedure that is 
standard in the study of deterministic bifurcations. The existence of a random contribution 
to both fast and slow time scales lies at the origin of the difference. 

We now extend the model presented in section |lj to include a random component in or. 
Physically, this corresponds to a random component in the control parameter of the system 
that has a significant frequency content at u <C ujh, and a correlation time that is small in the 
slow time scale emerging close to the bifurcation (inversely proportional to Or) . Under these 
conditions, it is also possible to assume that the random component is gaussian and white. 
For the three examples given in Section |, this stochastic component reflects, for example, 
the presence of fluctuations in either the temperature or gravitational field (Rayleigh-Benard 
or double-diffusive convection) or in the applied voltage (electro-hydrodynamic convection). 
Results pertaining to the onset of standing waves are presented in Section |111 A| while the 
transition to traveling waves is studied in Section [] 1B|. All numerical simulations reported 
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below have been performed using an explicit integration scheme, valid to first order in At 
(see Appendix). 



A. Bifurcation to standing waves 

Random fluctuations in the control parameter cir are introduced by letting a# — > + £ 
(with an now representing an average value). Since it is assumed gaussian and white, the 
noise £(i) obeys the statistics (£) = and (£(£)£(£')) = 2«5(t — £'), with k its intensity. In 
the region corresponding to the onset of standing waves (b > a.j), it is useful to introduce 
the variables A = x + y and Z = x — y, in terms of which Eqs. (j|) to (|j) become 

«9 4 A = (o B + bcos X )A + ^A(A 2 + Z 2 ) + 9 -fA(A 2 - Z 2 ) + £(t)A, (13) 

$Z = (o fl - bcos X )Z + fz(A 2 + Z 2 ) - 9 -fz(A 2 - Z 2 ) + £{t)Z, (14) 

and 

9 tX = 2az + ^(A 2 + Z 2 ) - 2&sin X ^^ 2 2 . (15) 

Linear stability analysis performed in Section |TT] showed that, in the deterministic case, 



standing waves appear supercritically at or = — d = —\Jb 2 — a 2 . From the linear part of 
Eq. (PH), this implies cosx = ^/l — a 2 /6 2 or x = X = arcsin(aj/6) at the (deterministic) 
bifurcation point. The trivial solution Z = remains stable above onset as the linear 



coefficient clr — bcosx in Eq. ([14] ) is negative. Thus, an initial difference between the 
amplitudes of the left and right traveling waves rapidly decays to zero. This is qualitatively 
unchanged in the stochastic equation, as the variable Z multiplies the noise and thus 
suppresses the influence of fluctuations as it goes to zero. Hence, close to onset, the governing 
equations can be approximated by 

d t A= (a R + bcos X )A+ 1 ^A 3 + £(t)A, (16) 

and 
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11; 



dtX = 2ai + — A 2 — 2b sin x- 



(17) 



Furthermore, in the weak noise limit, it is reasonable to expect the phase angle x to differ 
only slightly from its deterministic value at onset, % Therefore, we introduce the variable 
9 = x — X an d assume 9 << 1. Expansion of the trigonometric functions in equations ( |T6|) 
and flT7| ) yields to first order in 9 





A 




Or + d 




A 














+ 




9 




-2d 




9 





2f A 3 - A 



(18) 
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Just above onset, the two eigenvalues Ai = a^ + d and A2 
clr + d 

M = 

-2d 
unity except in the c 



-2d of the linearization matrix 
are of opposite sign. Furthermore, Ai << 1 while | A2 1 is of order 



ose vicinity of the codimension-two bifurcation point (where d — > 0). 
This implies the existence of two different time scales in the problem, the first one of which 
characterizes the rapid relaxation of the system to the center manifold 9 (A). To lowest 
order in A, 9 a (A) = {n i /Ad)A 2 ) as seen by letting d t 9 = in Eq. ([TSj) . The fact that £(t) 
does not appear in the equation for 9 allows us to neglect fluctuations away from 9 (A). 
The subsequent evolution of the system is therefore confined to the center manifold, and 
we look for stationary solutions of the form V(A,9) = P(A)5(9 — 9 (A)) to the Fokker- 
Planck equation corresponding to Eq. (|T8|). Explicitly, the time-independent probability 
distribution V(A, 9) describing the statistical properties of the system obeys the equation 



d_ 

dA 



(or + d)A + ^A 3 - ai9A + kA 



9 

v ' K dA 



A 2 V 



(19) 



d_ 

09 



-2,W + jA 2 + aj 2 ) V 



0. 



The fast variable is eliminated from the dynamics by integrating this equation over 9, with 
V(A, 9) = P(A)5(9 — 9 (A)). The second term on the L.H.S. vanishes once the integral is 
performed, as it is proportional to V(A, 9) evaluated at the limits of integration. This leaves 
an ordinary differential equation for P(A) 
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[a R + d)A + ( n R - ^-~kA 



P - kA 2 4tP = 0, (20) 
dA 



with solution 



P{A) =NA « — ^xp 



CLiTli \ j\. 



d J 8k 



(21) 



This distribution is normalizable (with M = 2[— (n^ — aiiii/d) /Sk] 2 ^ /T[(a,R + d)/2n] the 
normalization constant) as long as a R > —d. Below that value, P(A) = 6(A), which implies 
that, just as in the deterministic case, the value a R = —d marks the onset of standing waves. 
Just above onset, the expression given in Eq. (|2l|) exhibits a divergence at the origin (figure 
IDA). At clr = —d + k, this divergence transforms into a maximum which moves to the 
right as the control parameter is further increased (figure |2[B). Both figures, corresponding 
to a noise intensity k = 0.01, show excellent agreement between predictions from Eq. (|2T|) 



and the corresponding stationary distribution function obtained by integrating Eq. (18 
numerically. The simulations were performed using a time increment At = 0.01 and a 
bin size AA = 0.005. Initial conditions for A and 6 were chosen randomly from a uniform 
distribution in the interval [0, .05]. Results from 500 independent runs were used to compute 
P{A). Each run consisted of five million transient iterations after which a new point was 
added to the statistics every 500 iterations (for a total of 1000 points per run). 

From the distribution P(A), the various moments of A can also be determined. In 
particular, the standing wave's average amplitude is given by 



(A) 



ami \ 1 



n R 



/2 £(^+i) 

■p I a r +d \ 
\ 2k ) 



(22) 



d J 8k. 

As shown in figure QC, results from numerical integration of both the reduced set [Eq. 
(|18|)1 and the original equations for x, y and x are once again in excellent agreement with 
predictions from Eq. (|22|). As before, the values At = 0.01 and k = 0.01 were used in each 
of the 50 runs performed for each value of the control parameter. Each run consisted of 11 
million iterations (10 7 transient) with new points added to the statistics every 1000 time 
steps. At any value of the control parameter a R , (A) < Adet, with Adet the amplitude of the 
standing wave in the deterministic case (figure |||C) . 
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The statistics of the fast variable 9 also follow from the analysis given above. For instance, 
the average phase difference between the left and right components of the standing wave is 
given by 



r+oo r+oo r+oo 

(9) = / dA d9 9P(A)S(9 - 9 (A)) = / dA o (A)P(A) (23) 

JO J-oo JO 

(a R + d). 



4d 



ami \ 1 



1 '( ^ + 1) m 

a-H+ d \ d 
2k 



—, n R 

d 



Thus the average phase difference (9) or, equivalently, the average second moment (A 2 ), 
grows linearly with the control parameter a R . Furthermore, the slope characterizing this 
linear increase is independent of the noise intensity, so that both (9) and (A 2 ) assume their 
deterministic values. Figure |2[D compares predictions from Eq. ( |2"3"| ) with results from 
numerical simulations. 

The separation of time scales used to obtain Eqs. fl2~l|) and (^) gradually disappears as 
the co dimension- two point is approached from above (i.e., as | Ag | = d — *■ 0). Furthermore, 
since 9 D oc 1/d, fluctuations in the phase variable x grow in the vicinity of the TB point, 
implying that higher order terms in 9 should be kept in Eq. fll8|) . Although the predictions 
from Eqs. ([H]) and (22) fail when b ~ a^, numerical integration of Eqs. (^|) to (§) indicates 
no qualitative change in the system's behavior: the line marking the onset of standing waves 
from the trivial state remains un-shifted from its deterministic location, while the waves' 
average amplitude above onset is comparatively smaller. 



B. Bifurcation to traveling waves 

The transition to traveling waves can occur either from the trivial state (b < a,) or from 
a pre-existing standing wave pattern (b > a^). The TW state is characterized by a time- 
dependent phase angle (p — (fi+<f2 an d a finite difference in amplitude between the two wave 
components (Z ^ 0). Since A and Z evolve over similar time scales, the governing equations 
for x, y, and x can n °t be simplified as in Section |111 A|. However, results from a numerical 



study indicate that the general conclusions of Section |III A| also hold in the region below the 
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TB point. Hence, to the accuracy of the computations, no shift was detected in the location 
of onset. The latter was determined by computing the asymptotic amplitudes of the left 
and right traveling waves at different values of the control parameter or. For each one of 
these values, the complex equations (|3|) and (|J) (with noise included in or) were integrated 
numerically two billion times, using a time step of maximum size At = 0.0005. For all values 
of b < ctj considered, the bifurcation was observed at clr = ± 0.002, a value consistent with 



its deterministic location (or = 0). As in Section pll A| , a decrease in the traveling wave's 
average amplitude compared to the deterministic value was also noted above onset. 

The presence of fluctuations in the control parameter affects the emergence of TW above 
the TB point in a different way. The transition from SW to TW, which occurs along the 
oblique line in figure [I], in the deterministic case, takes place over a range of control parameter 
values when noise is added to the system. This smearing of the bifurcation is due to the 
fact that both states involved in the transition have associated amplitudes x and y which 
are non-zero. Therefore, contrary to the primary bifurcation, fluctuations contribute to the 
dynamics on both sides of the bifurcation point. Another way to see this is to define the 
variable x' = x — Xdet (and, similarly, y' = y — ydet), with Xdet the deterministic value of x 
at the bifurcation. The resulting equation for x' involves the stochastic term (xdet + %')£,{t), 
in which the noise multiplies both the small variable x' and the constant Xdet- The second 
component contributes additively to the dynamics, leading to an imperfect bifurcation. The 
transition from a standing to a traveling wave state therefore involves intermediate values of 
the control parameter a« for which the system behaves sometimes like a SW and sometimes 
like a TW. Numerically, this bifurcation interval was determined by monitoring the temporal 
evolution of the quantities Z and dt<fi which are both zero if the pattern is a SW. For the 
parameter values given above and for all the driving intensities considered, the interval was 
found to include the deterministic location of onset. 
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IV. STOCHASTIC VARIATION OF THE MODULATION AMPLITUDE 



If the random component in the externally controlled parameters has a significant fre- 
quency content around uh, the analysis given in Section [TT1] needs to be modified. We have 
first considered the case in which the external driving is a = (b + e 2%Uat . If the correla- 
tion time of £(t) is large compared with 1/ujh but short in the slow time scale emerging at 
the bifurcation, then £(£) can again be assumed to be gaussian and white. A more general 
choice of a would involve both a random amplitude and phase. In that case, the coupling 
coefficient in the normal form is no longer real and the analysis is somewhat more involved. 
The resulting stability diagram is qualitatively the same than the one presented below, and 
will be discussed elsewhere. 



A. Bifurcation to standing waves 



As in Section |III A| we first let 6 = 6 + ^ (t) and rewrite Eqs. @ to @ in terms of the 



variables A and Z. This gives 

d t A ={a R + b cos x)A + |i(i 2 + Z 2 ) + ^-A{A 2 - Z 2 ) + f (t) cos X A, (24) 

8 t Z = (a R -b cos X )Z + y Z ( A2 + ~ ^f z ( A2 ~ Z 2 ) ~ Z(t) cos X Z (25) 

and 

dtX = 2az + | (A 2 + Z 2 ) - 2&sin X ^±|J - 2£(f) sm X ^|J. (26) 

Close to onset (clr ~ —d, with b > aA, the variable Z quickly decays to zero and consequently 
drops out from the above equations. Thus, 

A 3 

d t A = (a R + b cos x)A + n R — + f (t) cos xA (27) 

and 

d tX = 2a, + -^A 2 - 2b sin X ~ 2£(t) sin x . (28) 
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Due to the presence of non-linear functions of A and x i n the stochastic terms, the Fokker- 
Planck equation corresponding to Eqs. (p7| ) and (pS| ) cannot be solved exactly. However, in 
the limit A — > 0, the term n{A 2 /2 on the R.H.S. of Eq. (^) can be neglected, effectively 
decoupling equation (^) from Eq. (^). Although this approximation is expected to hold 
only in a very small neighborhood around the bifurcation point, it is nevertheless sufficient to 
determine analytically the location of onset, which marks a transition from a state with A = 
to one in which (A) is arbitrarily small (although non-zero). The stationary probability 
distribution of the now independent variable \ obeys the following Fokker-Planck equation 

dP(x) ( (h b 



dx \ 2k sin 2 \ 2k sin x 

which yields 



-cot* }P(x) = 0, (29) 



[ tm(x /2)]-^ exp /_^y 
smx \2k tan x / 

This expression for P(x) is plotted in figure ^ for the average modulation amplitude b = 2.25. 
The distribution has a maximum close to x — arcsin(aj/6), a divergence at x = 71 an d a 
minimum at some intermediate value Xm (figure inset). Except when b w a iy P(Xm) « 1 
so that trajectories are most of the time confined to the interval [0,Xm]- Since the phase 
angle x evolves independently of A and over a much shorter time scale, it effectively acts 
in Eq. ( f27D as a second noise source, with non-zero correlation time and non-gaussian 



statistics. Rewriting the remaining equation for A using the variables (with zero mean) 
£' = cosx — (cosx) and £" = ^ cos% — cosx), we have 

A 3 

d t A = (a R + 6(cos X ) + (t(t) cos X ))A + n R — + b£'A + £"A, (31) 

which describes a pitchfork bifurcation taking place at 

[a R )thr = -b(cos x) ~ (f cos x) • (32) 

Using the Furutsu-Novikov theorem yi6| , the second average on the R.H.S. of Eq. 
simplifies to 
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(£cosx) = k(5cosx/SQ = 2K(sin 2 x), (33) 

so that [aijtfw = — b(cosx) — 2k (sin 2 x). Both averages are easily calculated from the 
probability distribution Eq. fl30|), normalized in the interval [0, Xm]- As shown in figure |], 
excellent agreement was found between predictions from Eq. (|32]) and numerical estimates 
obtained directly from Eqs. (Q), (0) and (|8|). In both cases, the location of onset is shifted, 
indicating a stabilization of the trivial state. Simulations were performed at the two values 
of cir delimiting each error bar in figure f|. The existence of a bifurcation within the interval 
was inferred from the large change in the asymptotic amplitudes and y^ noted across 
the interval. Ten runs were performed for each value of a# using a time step At = 0.005 and 
a total number of iterations per run N = 5 x 10 7 . Although figure ^ only shows results in 
the range 2.2 < b < 2.45, similar agreement was observed at larger values of the modulation 
amplitude. As mentioned above however, difficulties arise when b « (i.e., close to the 
TB point). To understand the origin of these difficulties, consider the temporal evolution 
of the phase angle x during a typical run at b = 2.25 (figure ^). Long periods during which 
X fluctuates according to the distribution Eq. fl30"| ) are followed by short intervals in which 
it rapidly increases by 2tt. The existence of such steps follows from the fact that P(Xm) is 
not identically 0, allowing trajectories in x-space to leave the interval [0, Xm] after a certain 
time. The average time a given trajectory takes to escape is given by the expression [IJ| 



1 f* m dv fy 
T = ^7 D , 2 / P(z)dz (34) 



2k J Xo P(y) sin y Jo 

which, evaluated for b = 2.25, yields T = 1.30 x 10 4 , a result essentially independent of 
the initial condition Xo- This estimate is in good agreement with the value T = 1.38 x 10 4 
obtained by averaging the number of jumps in x occurring during numerical simulations. 
In the case just considered, these jumps are rare and therefore statistically insignificant. 
Hence, for b = 2.25, Eq. ( |30| ) provides an accurate estimate for the averages (cosx) and 
(sin 2 x) which determine the location of onset. However, as the modulation amplitude b is 
lowered, P(Xm) increases and so does the number of steps in x- The analytical approach 
developed above eventually fails and the location of onset must be determined numerically. 



The black dots shown in figure 0, mark the onset of standing waves over the entire range 
of 6-values. They indicate a qualitative change in the system's behavior near the TB point 
with standing waves becoming stable with respect to the trivial state when b « dj. SW are 
also present in the region b < at, where they were previously unstable to traveling waves. 
Therefore, periods of rapid increase in x tend to favor the formation of a SW, the maxima 
and minima of which get inverted with each jump in x ( see Eq. fllPf) with x ~ * X + 27r). 
This stabilization of the SW state can be understood by first noting that close to the TB 
point, x = arcsinaj/6 ~ tt/2. Hence, in the deterministic limit, the terms proportional to 
cosx in Eqs.(||) and (0) go to zero, suppressing the mutual excitation between left and right 
traveling waves responsible for the emergence of a SW pattern. By increasing the probability 
of finding values of x away from x, the sudden jumps described above restore part of this 
constructive interaction and lead to the observed change in behavior. When the average 
modulation amplitude is small (b << 1), ~ 2a^, as seen by letting A, b — > in Eq. 
(p8|). From Eq. (|T0|), this result implies that, in the limit of small driving amplitude, the 
system is effectively oscillating at a new frequency w' H = wn + cii- Furthermore, the location 
of onset is at 

[aR\thr = -fe(cos(2at*)) - 2K(sin 2 (2ait)) = -k, (35) 

where the ensemble average has been replaced with an average over time. Results from 
simulations performed in the limit b << 1 agree with predictions from Eq. (|35"|). 



B. Bifurcation to traveling waves 

The empty circles in figure ([!]) mark the onset of traveling waves when b is a fluctuating 
quantity. As in Section [111 Bj these points were obtained from numerical simulations of Eqs. 



(0) and (H) in which changes in the amplitude difference Z and phase angle <fr were used to 
monitor the progressive transition from a standing to a traveling wave state. The results 
indicate a shift in the location of onset with the direction of this shift depending on the value 
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of the modulation amplitude. In particular, for values of b around or below the detuning dj, a 
delay in the onset of traveling waves was observed. Thus, the TB point, which is a distinctive 
feature of the deterministic stability diagram, disappears when a random component is added 
to the driving. 
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APPENDIX A: 



Numerical integration of the various stochastic differential equations encountered in sec- 



tions m and |V| was performed using an explicit scheme valid to first order in At. Expressed 
in terms of the Stratonovitch calculus, the algorithm used ||19[ maps the Langevin equa- 
tions 

Xi = fi({x k (t)}) + 9i({x k (t)})at), (Al) 
with £(t) gaussian white noise, to the discrete set 

Xi(t + At) = Xi (t) + fi({x k (t)})At + 9i ({x k (t)})E(t) (A2) 

The random number is gaussian distributed, with variance (S(t) 2 ) = 2nAt. As an 
example, we give the discretized version of Eqs.(f|) to (§) with noise included in the control 
parameter clr: 

x{t + At) = x{t) + At[(a R + bcosx(t))x(t) + c R x(t)(x(t) 2 + y{tf) + g R x{t)y{t) 2 } (A3) 
+x(t)E(t) + \x(t)~(t) 2 , 
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y(t + At) = y(t) + At[(a R + 6 cos X (t))y(t) + c R y(t)(x{t) 2 + y{tf) + g R y(t)x(t) 2 } (A4) 
+y(t)E{t) + ^(t)S(t) 2 , and 

X (t + At) = x (t) + At[2a, + n,(a;(t) 2 + y{tf) - bsm X (t)(x(t) 2 + y{t) 2 )/xy)\. (A5) 
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FIGURES 
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FIG. 1. The solid lines are the stability boundaries of Eqs. dH), (0) and (||) with 
Oj = 2,cr = — l,Cj = 2,gn = —1 and gi = 1. (•) : onset of standing waves when fluctua- 
tions of intensity k = 0.01 are added to the modulation amplitude b; (o) : onset of traveling 
waves under the same conditions. 
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FIG. 2. A and B: probability distribution for the standing wave's amplitude above onset for 
the same values of the parameters as in Fig 1, and, A, clr = —1.995; B, an = —1.98. C and D: 
average amplitude (A) and phase difference (9) as a function of the average control parameter clr. 
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X 

FIG. 3. Probability distribution V(x) corresponding to the average modulation b = 2.25. The 
distribution has a minimum at some value Xm above which it starts increasing again and eventually 
diverges at \ = ^ (inset). 
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FIG. 4. Location of onset in the presence of a fluctuating modulation amplitude b. For the 



average driving intensities shown, the bifurcation point is shifted to the right of its deterministic 



position. 
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FIG. 5. Temporal evolution of the phase difference x during a typical run at b = 2.25. 
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